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Abstract 

We start with the QCD sum rules which are originally based on the idea that it is the power-like corrections to the 
parton model which are related to the confinement. The naive use of the Operator Product Expansion ensures that 
there is a 'gap' in the powers of Aqcd which miss the quadratic terms and start with the quartic term, proportional 
to the gluon condensate, < (G^ v ) 2 >. We review how this hypothesis stood against various checks through the last 
three decades and how it was modified through inclusion of the missing link, that is quadratic corrections. In field 
theoretic language the quadratic corrections are dual to long perturbative series. In the dual description, the quadratic 
corrections are conveniently parameterized in terms of the metric in extra dimensions. We emphasize that the dual 
models do not incorporate the so called infrared renormalon. 
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1. Introduction 

This talk was given at a special session of the QCD 10 
conference at Montpellier, devoted to (30+1) year of the 
QCD sum rules [1]. I am thankful to Stephan Nari- 
son for organizing this event and, also, for overtaking 
the most difficult part of the job, that is reviewing the 
present status of the sum rules |2|]. As for myself, I 
chose for myself to talk more about not yet fully settled 
issues. The emphasize is mainly on the quadratic cor- 
rections, simply missing from the standard sum rules. 
To a large extent the talk is based on the published pa- 
pers JslHl but we add a few remarks as well. 

Originally the sum rules were applied to two-point 
functions induced by external currents. More specifi- 
cally, one considers integrals of the kind 

f(M 2 ) = j dsexp(-s/M 2 )ImTl(s) , (1) 

here M 2 is a large mass parameter, /miles') is imaginary 
part of a polarization operator and, speaking generically, 
can be measured by studying transition induced by cur- 
rents. The sum rules equate the observable (Q]i to power 
corrections: 

A OCD 

f{M 2 ) « [f(M 2 )] parton model(l + c gl.cond—j^-) (2) 
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where the product c s ;. con dAi C£) is calculable in terms of 
the gluon condensate, or matrix element < (G^ v ) 2 >. 

The form (f2]) is oversimplified since it omits, in 
particular, the first perturbative correction of order 
a s (M 2 ) ~ In -1 M 2 , the quark-condensate terms and so 
on. But nevertheless, eq. (f2]l does summarize correctly 
the basic idea that in the crucial region of M 2 it is the 
power corrections which are related to the confinement 
(manifested through resonance masses and widths en- 
tering (Q}). The power corrections are given in terms 
of matrix elements of local gauge invariant operators. 
In particular, there is no correction of order Ag CD since 
there is no corresponding operator of dimension d=2. 
Eq. (f2]i is in no way obvious and was introduced on 
phenomenological grounds Jll . In particular, dropping 
the perturbative corrections which are powers of a s (M 2 ) 
and keeping the power correction which is of order 
exp(-const/a s ) might look confusing. 

2. Quadratic corrections: 

unifying continuum- and lattice-languages 

The expansion (f2]l works well in many cases J2l- 
However, there exist reasons to revisit it: 

1 . Trying to improve the quality of the sum rules one 
starts to calculate perturbative corrections. The 
question is, whether the numerical value of < G 2 > 
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is kept independent of these corrections or should 
it vary from order to order of perturbation theory? 
2. The Cornell potential for heavy quark interaction: 



V QQ (R) 



const 



R 



+ crR , 



(3) 



where R is the distance between the quarks and cr is 
the string tension holds numerically at all the dis- 
tances measured on the lattice including small dis- 
tances. In fact, this is an ideal example of a kind 
of expansion ©, with a single power-like correc- 
tion present. The problem is that the correction to 
the leading, Coulomb-like term at short distances 
is quadratic, ~ (crR 2 ). And quadratic corrections 
are not originally included into © and things look 
mysterious. 

Although the questions are straightforward, convinc- 
ing answers are difficult to get phenomenologically. The 
reason is that, after all, we talk about relatively small 
power-like corrections. Each particular case of such 
phenomenology is tedious, and more important, very 
difficult to follow from outside. Nevertheless, after a 
few years of hard work the answer seems to be unique: 
quadratic corrections do exist • 

Moreover, the quadratic correction is directly related 
to confinement. Although the statement might look too 
strong and vague, actually, it has a well defined con- 
tent. Namely, on the lattice one is able to clarify what 
kind of field configurations are responsible for the con- 
finement. In the lattice nomenclature these configura- 
tions are called monopoles and vortices, for a review 
see, e.g., |g|. Most amusing, they occupy a small frac- 
tion of the lattice which tends to zero with the vanish- 
ing lattice spacing a — > 0. Closer to our story, it was 
demonstrated that the non-perturbative, i.e. monopole- 
or vortex- related potential is indeed linear at all the dis- 
tances beginning with a single lattice spacing: 



\ in 



cr R 



(4) 



/ non-perturbative 

Details and references can be found in 0], Thus, the 
short-distance quadratic correction is certainly there ! 

3. Power corrections vs perturbative series 

What is the relation between the power corrections 
and perturbative series? Concentrate on the gluon con- 
densate itself. Then perturbatively: 

° 4 T^(^)<7 GG ^ S ".( 1+ °"<)' (5 » 



where a s = g 2 /4-n is the strong interaction coupling, 
/3(g) is the beta function, a n are perturbative coefficients. 
Moreover, a is the lattice spacing and the factor a 4 in 
the l.h.s. of Eq © is introduced to cancel the UV diver- 
gence inherent to the quantity considered. 

The perturbative expansion (O is expected to be 
asymptotic. Namely, for n large enough the expansion 
coefficients are expected to grow factorially: 



lim a„ 

n—>oo 



- »'(t)" 



(6) 



The series (|6]l is called infrared renormalon. If we es- 
timate the uncertainty of the asymptotic expansion due 
to the factorial growth (O we find contribution of order 
(Aqcd ' fl ) 4 compared to the leading perturbative term. 
Such an uncertainty would introduce < G 2 >~ Ag CD 
which is the physical gluon condensate entering, in par- 
ticular QCD sum rules (f2j). This is the standard wisdom 
on the relation between divergences of the perturbative 
series in large orders n and the Operator Product Expan- 
sion. 

Where is then the hypothetical quadratic correction? 
Well, there is no pronounced role for such a correction 
in the set up considered. It should be buried within 
the still-convergent orders of the perturbative series. In 
other words, keeping the quadratic correction might be 
reasonable only as far as the perturbative series is not 
long enough. If we keep many terms, then the quadratic 
correction is to be eaten up by the perturbative terms. 
Note that within such a logic the quadratic correction 
(if any) is inferior to many orders of perturbation se- 
ries and uninteresting. Moreover, the very language of 
power corrections seems rather irrelevant if we need to 
keep explicit many orders of perturbation theory. 

Things become, however, much more interesting if 
we turn to the example of the longest perturbative series 
known. It is indeed for the gluon condensate <j5j and 
contains 20 (no mistake: twenty) first terms in the ex- 
pansion, see |Ht] and references therein. Moreover, the 
full value of the gluon condensate is known from the 
lattice measurements since it is simply the plaquette ac- 
tion and can be measured to a very high precision. As a 
result one can use the following fitting procedure: 



(AP) N = P full - Z n p„a" s » (A QCD ■ a) p(A ° 



(7) 



where P full is the exact (or full) plaquette action, p„ are 
perturbative coefficients evaluated explicitly and the dif- 
ference between the partial (up to order N) perturbative 
series and the full value is fitted by a power-like correc- 
tion where the index p(N) depends on N itself. 



The results [8] concerning the index p(N) are remark- 
able. Namely, the fits produce: 



p(N) « 2 i/JV < 10 
p(N) ~ 4 i f N > 10 



(8) 



Thus the perturbative series does know about the sacred 
quadratic and quartic corrections although this knowl- 
edge is very difficult to express analytically. Everything 
is numerical. 

Another remarkable result is the simplicity of the co- 
efficients p„. The expansion is close to a geometric se- 
ries with the ratio 



Pn+l 
Pn 



Const . 



(9) 



and we do not quote the numericalvalues of the fit pa- 
rameters here. Observation (01 implies that there is no 
sign of the infrared renormalon at least in the first 20 
orders of expansion. This is the reality which we have 
to confront and appreciate theoretically. 

In conclusion of this section it is worth mentioning 
that similar interplay between perturbative and power- 
like corrections was observed earlier by analyzing much 
shorter perturbative expansions. For example, it was ob- 
served in Ref. |9[] that the quadratic correction at short 
distances to the Coulomb-like potential (see Eq. ©) is 
reproduced in fact by higher orders of perturbation the- 
ory. Another example is provided by Ref. [10] where 
sum rules for the structure function xF?,(x) measurable 
in deep inelastic neutrino scattering were analyzed. It 
turned out that including higher orders of the perturba- 
tive expansion reduces considerably the numerical value 
of the quadratic correction. Note that in case of the 
xFi(x) function the quadratic correction is commonly 
associated with the infrared renormalon while in case of 
the quark potential the quadratic correction corresponds 
to the ultraviolet renormalon Q. Strong correlation be- 
tween the quadratic correction and perturbative series is 
found in both cases. The limitation of the analyses just 
mentioned is that the perturbative series known explic- 
itly are relatively short, say 3-4 terms, and to reach con- 
clusions one relies on theoretical estimates of higher- 
order terms. Signals for importance of quadratic correc- 
tions were obtained in quite a few other analyses, see in 
particular 111 lh . 



1 Mostly, when we talk about the quadratic corrections we have in 
mind vacuum-to-vacuum matrix element of a correlator of two cur- 
rents. Then the quadratic correction might come only from short dis- 
tances and has no straightforward interpretation in terms of the OPE. 
In case of the deep inelastic scattering, as is well known, quadratic 
correction is associated with large distances. 



To summarize, there is strong evidence in favor of 
crucial role of the unconventional quadratic corrections 
in correlators of two currents. Moreover, the quadratic 
correction is dual to a long perturbative series. Thus, 
one is to use either of them but not both. "Unconven- 
tional" means that there is no operator of the dimension 
d = 2 which would enter the OPE. Thus, it is difficult 
to even parameterize the quadratic correction within the 
field-theoretic approach: the quadratic correction is a 
part of the coefficient function in front of the unit oper- 
ator but we are lacking means to distinguish it from the 
rest of the coefficient function. 

4. Dual models and power-like corrections 

4.1. Finding a dimension-two parameter 

The situation with the quadratic correction is quite 
paradoxical. On one hand, there is accumulating and 
strong phenomenological evidence in favor of such a 
correction. On the other hand, the field-theoretic frame- 
work does not even provide us with a suitable parame- 
terization of the quadratic correction. The most success- 
ful phenomenological model [5] introduces a "short- 
distance gluon mass" m g , so that one replaces the gluon 
propagator by 

q L q L + m~ 

where the gluon mass turns to be tachyonic. Clearly, the 
replacement ( flOb is consistent with the gauge invariance 
only in the Born approximation. In higher orders, there 
is no way to introduce a dimension-two quantity. 

The resolution of the paradox seems to be that the 
quadratic correction belongs rather to the realm of the 
dual models for the Yang-Mills theories. As is well 
known, the dual or stringy formulations of non-Abelian 
gauge theories is derived only in some supersymmet- 
ric cases. In QCD case, there are educated guesses 
on the universality class 11121 1 1311 and more explicit but 
rather arbitrary models, see, in particular, B14I416I1 . The 
main point now is that the models allow to introduce a 
dimension-two parameter to the theory in a gauge in- 
variant and natural way. 

To begin with, the stringy models readily incorporate 
confinement. Indeed, linear potential at large distances, 
see is commonly related to existence of a string con- 
necting the quarks. To be more ambitious and describe 
the potential at all the distances one can introduce a run- 
ning string tension cr(l) where I is the length of the string 
and cr{l) is its tension as function of the length. More- 
over, one can trade the running tension for a string liv- 
ing in an extra dimension z with non-trivial geometry, 
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G mn (z). The action of the string is the Goto-Nambu ac- 
tion, 

S = 2^7 J d 2 ^detG nm d a X^d p X n . (11) 

The relation to the potential is provided by the boundary 
condition that the string ends on the Wilson line 1 1 811 . 



< W(C) > = exp(-S(C) mi „), 



(12) 



and the Wilson line belongs to z — 0. One can then 
trade the potential (]4]l for an explicit form of the metric 
G m „(z). In particular, the following metric reproduces 
the Cornell potential (|4]i quite satisfactory II 1611 : 

cz 2 /2 

ds 2 = G mn dX m dX n = R 2 ——(dx'dx' + dz 2 ), (13) 

z z 

where the ordinary 4 dimensional space corresponds to 
z = 0, R 2 is a constant, the parameter c is of phenomeno- 
logical nature and its numerical value, c « 0.9 GeV 2 , is 
obtained from the fitting procedure. 

Most remarkably, the parameter c specifies the 
quadratic correction at small z in a perfectly gauge in- 
variant way. Thus, the metric ( TT31 allows to introduce 
the quadratic correction in the universal geometric lan- 
guage. Upon fixation of the coefficient c in terms of 
the potential one can, in principle, evaluate quadratic 
corrections to other observables. For technical reasons, 
there are not many examples of this kind. However, one 
can say that the emerging phenomenology turns to be 
successful, see, in particular 1 15 , 1^1 . 

To summarize, the use of the dual models promotes 
the numerical observation on the dominance of the 
quadratic correction at small R in the heavy-quark po- 
tential (0]l to the status of a universal feature of the short- 
distance physics. 

4.2. On the sign of the quadratic correction 

There is a puzzling question about the sign of the co- 
efficient c in the metric ( fob . We fixed it as positive, by 
mapping the potential (0]l into the metric. Actually, the 
sign is fixed merely by the condition that we should have 
confinement, as first observed in lfl6Tl . see also lfl7ll . On 
the other hand, papers which introduced the quadratic 
correction originally Il4l 1 1 511 are using rather negative 
value of the coefficient c, h ~ exp(-cz 2 /2), instead of 
h ~ exp(+cz 2 12) as in Eq. ( TT31 ). 

This discrepance might be finally settled in favor of, 
say, positive sign of c through further phenomenologi- 
cal analysis. One cannot rule out, however, that there is 
a matter of principle behind this apparent discrepancy. 
Namely, the papers lfl4l [l5ll refer rather to Minkowski 



space while numerically impressive successful applica- 
tions MM 

of the metric refer to the Euclidean space. 
Then we would come to an intriguing possibility that in 
theories with confinement there is no analytical contin- 
uation of the quadratic correction from the Minkowski 
to Euclidean space. 

The question is then, whether the insight we got on 
the nature of the quadratic correction as being dual to 
long perturbative series supports the idea of the non- 
analyticity or rather rejects it. Analytic properties of the 
running coupling are discussed in many papers, see, in 
particular, [20] and it seems obvious that the continua- 
tion of the coupling from the Minkowski to Euclidean 
space and vice verse is not straightforward at all. 

To put it even simpler, by continuing the log factor in 
the running coupling one gets an imaginary part in the 
Minkowski region, In Q 2 = \Ti(-Q 2 ) + in. The imaginary 
part is big numerically and summing up higher orders 
in perturbation series at Q 2 > and at q 2 = -Q 2 might 
well result in different quadratic corrections. 

Basing on the example of the perturbative evaluation 
of the gluon condensate discussed above it is tempting 
to assume that geometric-series-type of expansion does 
incorporate the quadratic correction in other cases as 
well. And, indeed, one observes [3] that existing exam- 
ples of relatively long explicit perturbative expansions 
do not rule out that geometric series approximately ap- 
ply in other cases as well. For example, using the results 
of Ref. l2~lll one finds for the polarization operator as- 
sociated with the scalar currents in QCD: 



- Q2 dW 2mQ2) = ^( 1+5 - 67a 



d 

+45.85^ + 465.8a*;) + 5589arjj . 



(14) 



We see that this expansion numerically is quite close to 
a geometric series. In the context of our discussion, it 
is crucial that the expansion (TBI refers to the Euclidean 
coupling constant a s (Q 2 ). If one would use a s (-Q 2 ) 
as an expansion parameter, no similarity to a geometric 
series would persist. 

To summarize, the duality between the power cor- 
rections and long perturbative series rather supports the 
suggestion that the quadratic correction cannot be triv- 
ially continued from the Minkowskian to Euclidean do- 
main and vice verse. Much more is to be done, however, 
to strengthen the argumentation. 

4.3. Absence of infrared renormalon 

Probably one of the most important changes brought 
by the dual models to the dogma of high-order correc- 
tions is the absence of the infrared renormalon 0]. 
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In more detail, to imitate the gluon condensate (0 
in the continuum theory one evaluates vacuum expec- 
tation value of a circular Wilson line of small radius r, 
r 2 Ag CD <k 1 . In the dual model one uses eqs $12\ . (fTTl) . 
(TT3b to perform the explicit calculation. The result can 
be represented as: 



< W(r) >x; exp 



consul -0.84cr z 



-0.035cV 



■(15) 



The first two terms in the r.h.s. of (TT3T > correspond to 
the quartic and quadratic ultraviolet divergences in (0. 
The evaluation of these terms in the continuum theory 
is not unique and depends in fact on the regularization 
procedure. These terms cannot be compared directly 
with the results of the perturbative calculation (0 in the 
lattice regularization. 

The most interesting result is the term proportional to 
c 2 in the r.h.s. of Eq dlS) . It corresponds, in the contin- 
uum language to < G 2 >~ Ag CD and is most relevant to 
the sum rules (0 and phenomenological fits. Moreover, 
it is this term which is commonly related to the expected 
IR-renormalon divergence of the perturbative series (0. 
The crucial observation |4j is that in the dual-model ap- 
proach the < G 2 >~ Ag CD term is related in fact to short 
distances and is calculable. It corresponds to the z 4 term 
in the expansion of the factor h(z) = exp(cz 2 /2) in the 
metric dT3b at small z, z — » 0. 

This is quite a remarkable possibility revealed by the 
dual models: the < G 2 >~ Ag CD piece in the gluon 
condensate is associated with the small, not large dis- 
tances. If the piece < G 2 >~ Ag CO is indeed associated 
with short distances then it should be calculable in per- 
turbation theory. This prediction of the dual models is 
in perfect agreement with the absence of the infrared 
renormalon from the explicit calculation (0, see for 
discussion above. The ifference between the standard 
field theoretic estimates which indicate the presence of 
the IR renormalon and dual approach which, does not 
incorporate the infrared-sensitive piece of < G 2 > is 
of pure geometric nature. In the both approaches, the 
total value of the small Wilson line is dominated by 
small distances, of order r. The infrared sensitive piece 
is anyhow relatively small and is due to rather excep- 
tional or suppressed configurations which involve dis- 
tances much larger than r. One can readily check that 
it is much easier to reach such distances for one virtual 
particle (gluon) than for a string, see (TT2l >. ( TT3T >. 

4.4. Choice of the dual model 

We can summarize our discussion by saying that the 
mystery of the quadratic correction is resolved only by 
the dual models. Moreover, one can reverse the logic 



and try to clarify which dual models are fitting the 
knowledge on the power corrections in the best way. 
The conclusions then are: 

1 . the quadratic correction is to be built explicitly into 
the dual model 

2. the model is to be formulated in terms of strings, 
not just fields living in the extra dimensions 

It is interesting to note that an independent and much 
more thorough analysis of phenomenological implica- 
tions of the existing dual models results in a similar con- 
clusion Il22ll . Unfortunately, stringy models are more 
difficult to apply to phenomenology of the power cor- 
rections and mostly one considers models with fields 
(not strings) living in extra dimensions, see, e.g., 112 311 . 



5. Conclusions 

To summarize, there is substantial progress in under- 
standing of the power corrections brought to light by 
the sum rules. Namely, the emphasize shifted to the 
quadratic correction absent from the original, simpli- 
fied form of the sum rules. It turned out that it is just 
this correction which is most closely related to the con- 
finement. Moreover, introduction of this correction al- 
lows for a straightforward interpretation of some lattice 
data and unifies the continuum-theory and lattice lan- 
guages. Thus, phenomenologically the quadratic cor- 
rection resolves quite a few puzzles. However, interpre- 
tation in the field-theoretic language is rather awkward: 
it is dual to a long perturbative series, as is confirmed by 
a number of perturbative calculations, see in particular 
|f8|-|Toh . With the advent of the dual models of QCD the 
quadratic correction found its interpretation in terms of 
the metric in an extra coordinate z. Dual models which 
incorporate this quadratic correction at small z turn to 
be successful phenomenologically, probably even most 
successful. 
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